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Abstract 

A simple model based on relativistic geometry and final-state hadronic rescattering is used to 
predict pion source parameters extracted in two-pion correlation studies of proton-proton collisions 
at y/s = 14 TeV. By comparing the results of these model studies with data, it might be possible 
to obtain information on the hadronization time in these collisions. As a test of this model, 
comparisons between existing two-pion correlation data at y/s = 1.8 TeV and results from the 
model are made. It is found at this lower energy that using a short hadronization time in the 
model best describes the trends of the data. 
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I. INTRODUCTION 



With first proton-proton collisions at y/s = 14 TeV from the Large Hadron Collider (LHC) 
being only about a year (or so) away, it is tempting to use simple models to make baseline 
predictions of what we might expect for "bread and butter" observables at this unexplored 
energy. Comparisons between data and such models could give us a first impression of 
the presence of new physics which might cause significant disagreements between them. 
If significant disagreements are seen, the simple models might then be used to point in a 
direction as to the nature of the new physics. 

The "bread and butter" observable studied in the present work is two-pion correlations. 
From this observable, information about the space-time geometry of the pion emissions 
produced in the proton-proton collisions can be, at least in principle, extracted using the 
interferometric technique pioneered by Hanbury Brown and Twiss (HBT) [lj and first used in 
particle physics mp — p collisions by Goldhaber, Goldhaber, Lee and Pais (GGLP) jj]. Many 
such experimental studies using this technique have been carried out over the past nearly 50 
years 3|, |4| , the highest energy study being carried out at the Tevatron with p — p collisions 
at yfs = 1.8 TeV [5J. The strategy of the present study is to use a simple model based on 
relativistic geometry and final-state hadronic rescattering to predict pion source parameters 
extracted in two-pion correlation studies of proton-proton collisions at y/s = 14 TeV. As a 
test of this model, comparisons with existing two-pion correlation data at \fs = 1.8 TeV 
are made. A similar study which served as the inspiration for the present work has been 
published by Paic and Skowronski [6]. The main differences between that study and the 
present approach are that in the present approach 1) a somewhat simpler geometric picture 
of hadronization is used, e.g. no explicit identification of jets vs. non-jets is made, 2) 
for simplicity only 1-dimensional invariant correlation functions are studied, 3) at Tevatron 
energy Gaussian fits to the model-generated two-pion correlation functions are made to 
directly compare with experiment whereas at LHC energy the two-pion correlation function 
is fit to a more general function , and 4) the effects of final-state hadronic rescattering are 
included since particle multiplicities become relatively large at these higher energies. 

The paper is divided into the following sections: Section II gives a description of the 
model, Section III presents results of the model and discussions for p—p collisions at y/s = 14 
TeV and p — p at 1.8 TeV and comparisons with Tevatron data, and Section IV gives 
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conclusions and summary. 



II. DESCRIPTION OF THE MODEL 



The model calculations are carried out in four main steps: A) generate hadrons in p — p 
and p — p collisions from PYTHIA 3] , B) employ a simple space-time geometry picture 
for hadronization of the PYTHIA-generated hadrons, C) calculate the effects of final-state 
rescattering among the hadrons, and D) impose Bose-Einstein correlations pairwise on pions, 
calculate correlation functions, and fit the correlation functions with Gaussian or more 
general functions to extract pion source parameters. These steps will now be discussed in 
more detail. 



A. Generation of the p — p collisions with PYTHIA 

The p— p and p—p collisions were modeled with the PYTHIA code version 6.326. The 
parton distribution functions used were the same as used in Ref. jsj]. Events were generated 
in "minimum bias" mode, i.e. setting the low-pr cutoff for parton-parton collisions to zero 
(or in terms of the actual PYTHIA parameter, ckin(3) = 0). Runs were made both with 
y/s = 1.8 and 14 TeV to simulate Tevatron and LHC (full energy) collisions, respectively. 
Information saved from a PYTHIA run for use in the next step of the procedure were the 
momenta and identities of the "direct" (i.e. redundancies removed) hadrons (all charge 
states) it, K, p, n, A, p, u, r], rf ', 0, and K*. These particles were chosen since they are the 
most common hadrons produced and thus should have the biggest effect on the two-pion 
correlation functions extracted in these calculations. 



B. The space-time geometry picture for hadronization 

The simple space-time geometry picture for hadronization consists of the emission of a 
PYTHIA particle from a thin uniform disk of radius 1 fm in the plane transverse (x-y) to 
the beam direction (z) followed by its hadronization which occurs in the proper time of the 
particle, r. The space-time coordinates at hadronization in the lab frame (xh,Vh, Zh,th) for 
a particle with momentum coordinates (p x ,P y ,Pz), energy E, rest mass m , and transverse 
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disk coordinates (xo,yo) can then be written as 



x h = x + T— (1) 
m 

Vh = yo + t^- (2) 
m 

z h = r— (3) 
m 

E , , s 

t h = t— 4 
m 

The simplicity of this geometric picture is now clear: it is just an expression of causality 
with the assumption that all particles hadronize with the same proper time, r. A similar 
hadronization picture (with an initial point source) has been applied to e + — e~ collisions 9] . 
We do not a priori know the value of r but from the geometric scale of a p — p collision we 
might guess that r falls in the range < r <~ 1 fm/c. In order to study the dependence 
of the results of the model on this parameter, calculations will be carried out with a range 
of values. Note that the HBT results given later from the model are found to not strongly 
depend on the choice of the radius of the initial transverse disk within a range of 1 ± 0.5 fm 
or on the choice of a disk versus a smoothly dropping off distribution such as a Gaussian 
due to the effects of these assumptions being "washed out" by the randomizing effects of 
the "causality term" in Eqs.(l) and (2) and of final-state rescattering. 

C. Final-state hadronic rescattering 

Since very high energy p — p collisions are being considered here and the most interesting 
collisions are normally those producing the highest particle multiplicities, it seems possible 
that at early times during the collision the particle density could reach a level at which 
significant final-state hadronic rescattering might take place. An attempt is made to take 
this effect into account in the present calculations. 

The hadronic rescattering calculational method used is similar to that employed in previ- 
ous calculations for heavy-ion collisions at CERN Super Proton Synchrotron (SPS) energies 
and BNL Relativistic Heavy Ion Collider (RHIC) energies J, 111]], where particle densities 
are high enough to produce significant rescattering effects. Rescattering is simulated with 
a semi-classical Monte Carlo calculation which assumes strong binary collisions between 
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hadrons. Relativistic kinematics is used throughout. The hadrons input into the calculation 
from PYTHIA are pions, kaons, nucleons and lambdas (tt, K, N, and A), and the p, u, r], 
■)]', 0, A, and K* resonances. For simplicity, the calculation is isospin averaged (e.g. no 
distinction is made among a 7r + , 7t°, and 7r _ ). 

The rescattering calculation finishes with the freeze out and decay of all particles. Starting 
from the initial stage (t = fm/c), the positions of all particles in each event are allowed 
to evolve in time in small time steps (At = 0.1 fm/c) according to their initial momenta. 
At each time step each particle is checked to see a) if it has hadronized (t > t^, where th is 
given in Equation 4.), b) if it decays, and c) if it is sufficiently close to another particle to 
scatter with it. Isospin-averaged s-wave and p-wave cross sections for meson scattering are 
obtained from Prakash et al. 12| and other cross sections are estimated from fits to hadron 
scattering data in the Review of Particle Physics[13]. Both elastic and inelastic collisions 
are included. The calculation is carried out to 20 fm/c which allows enough time for the 
rescattering to finish (as a test, calculations were also carried out to 40 fm/c with no changes 
in the results). Note that when this cutoff time is reached, all un-decayed resonances are 



allowed to decay with their natural lifetimes and their projected decay positions anc 



times 



are recorded. The rescattering calculation is described in more detail elsewhere [4, |ll|. The 
validity of the numerical methods used in the rescattering code have recently been studied 
and verified Jjj] . 



D. Correlation function calculation and fitting 

For the two-pion correlation calculations, the two-pion correlation function is formed and 
either a Gaussian or more general function is fitted to it to extract the final fit parameters. 
In the present calculation boson statistics are introduced after rescattering using a method 



of pair- wise symmetrization of bosons in a plane- wave approximation [IJjj]. The final step 
in the calculation is extracting fit parameters by fitting a parameterization to the Monte- 
Carlo-produced two-pion invariant correlation function, C(Qi nv ), where Qi nv is the invariant 
momentum difference defined as the magnitude of the difference between the four-momenta 
of the two pions, i.e. Qi nv =| pi — P2 |. The forms of the Gaussian and general fit functions 
are given, respectively, by 
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C{Q mv )=A[l + \eM-R 2 QU] (5) 

or, 

C(Q mv ) = A[l + A cos (aQ 2 inv ) exp(-R 5 Q 5 mv )} (6) 

where R is a radius parameter, A is an empirical parameter normally employed to help 
fit the function to the correlation function (i.e. A = 1 in the ideal case of pure Bose-Einstein 
correlations), a describes oscillations in the correlation function, 6 represents the degree to 
which the correlation function falls off with increasing Qi nv , and A is a normalization factor. 
For the Gaussian simple connection can be made between R and the space-time 

distribution of the pion source, p(r), where r is a position variable, via 

(7) 

and, 

C{Q mv ) ~ 1 + A | ~p{Q inv ) | 2 (8) 

where p(Qi nv ) is the Fourier Transform of p{r) in terms of Qi nv . Inserting the Fourier 
Transform of Eq.(7) into Eq.(8) gives Eq.(5). The Gaussian function was used by E735 to 
extract R and A from data and thus is used exclusively to extract these parameters in the 
present calculations for p — p at t/s = 1.8 TeV to compare with the data. The general 
fit function, Eq.(6), is used exclusively to extract fit parameters to the model correlation 
functions for calculations of -y/i = 14 TeV p — p collisions. This is done since there are 
no restrictions on the fit function which can be used, the goal being to characterize the 
correlation function with as good a fit as possible to study the dependencies of the fit 
parameters on various kinematical conditions. The motivation for the particular form of Eq. 
(6) is discussed in detail elsewhere^]. Note also that a form similar to this has been used 



to fit preliminary pion correlation functions obtained in 



;he LEP L3 e + — e experiment in 



which a hint of a baseline oscillation has been observed 10]. 



III. RESULTS AND DISCUSSION 



Results from the model calculations described above are now presented and discussed. A 
comparison between model calculations for p—p collisions at = 1.8 TeV and experimental 
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results from the Tevatron E735 experiment are presented first as a reality check on the 
present model near the highest energy collisions presently available, followed by predictions 
from the model for p — p collisions at \fs = 14 TeV. 



A. Comparisons with data for p — p collisions at y/s = 1.8 TeV 

Although experimental two-pion HBT results for p — p collisions at \fs = 14 TeV are not 
yet available from the LHC with which to compare the predictions which will be presented 
later in this work, it is possible to compare results calculated from the present model with 
existing experimental two-pion HBT results from Tevatron experiment E735[5], which stud- 
ied p — p collisions at \fs = 1.8 TeV. Such a comparison with p — p collision data near the 
highest existing energy will point towards what expectations we should have for the present 
simple model to predict the higher LHC full-energy HBT behavior. 

To carry out this comparison, calculations were made with the present model for \fs = 1.8 
TeV p — p collisions using the same parton distribution functions in PYTHIA as for the 
yfs = 14 TeV case as mentioned above (this was done to be as consistent as possible 
with the \fs = 14 TeV calculations - it is not expected that using different pdf's in the 
model calculations would effect the present results significantly) . Gaussian fit parameters 
were extracted from the calculations using Eq. (5) since this was essentially the same fitting 
procedure used by E735 to extract the fit parameters R and A. The E735 parameters 
with which comparison is made in the present work were obtained directly from Table II of 
Ref. [5I for the Nq (see below) dependence and from Table III in the same reference (using 
their conversion R = 0.254 + 1.023-Rg, where Rq is defined in the reference) for the pt 
dependence. The E735 pion acceptance was simulated in the model calculations with simple 
kinematical cuts on rapidity and px- Dependency on the charged particle multiplicity in the 
E735 hodoscope, Nc, was also studied, being simulated in the model calculations with an 
acceptance cut[5J]. 

Figures [1] - [8] present results of the model for \fs = 1.8 TeV p — p collisions. Comparisons 
with E735 Gaussian fit parameters are shown in Figures H] - [HJ 

Figure [T] shows model rapidity and pr distribution plots for all final-state particles (i.e. 
pions, kaons, and nucleons) for three cases: 1) directly from PYTHIA, 2) PYTHIA with 
a Tj — pt "hole", and 3) same as 2) but with rescattering turned on and r = 0.1 fm/c. 
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For case 1), PYTHIA events are directly run through the model code without any other 
process applied to them except to decay As, and the p, u, rj, rf, 0, A, and K* resonances, 
i.e. "pure PYTHIA". Since PYTHIA has been tuned to agree reasonably well with exist- 
ing experimental data, including Tevatron data, these distributions should remain at least 
approximately the same after rescattering has been turned on. This turns out not to be 
the case for r = 0.1 fm/c, since It is found that if "pure" PYTHIA events are input into 
the calculation with rescattering turned on, a small peak results in the rj distribution near 
midrapidity and the pt distribution is overly enhanced compared with "pure" PYTHIA. 
This gives the first suggestion that final-state hadronic rescattering can play a noticeable 
role in these collisions. In an effort to compensate for the rescattering effects so as to give 
approximate agreement with the "pure" PYTHIA distributions, a rj — pt "hole" is inserted 
in the input PYTHIA events before rescattering, as shown, and the rescattering then fills 
the "hole" to approximately agree with the "pure" PYTHIA distributions, also shown in 
Figure [H For this case, the "hole" is defined by randomly throwing away 5% of the particles 
in the input PYTHIA events in the y — px region — 1 < y < 1 and px > 0.5 GeV/c. For 
the larger values of r studied, i.e. r = 0.5 and 1.0 fm/c, less rescattering takes place due to 
the larger initial hadronization volume and thus lower initial particle density and the "hole" 
depth is reduced, using the prescription 5%(0.1/r) to define it. The justification for using 
this "hole" method is that reasonable agreement with the "pure" PYTHIA distributions is 
obtained with recattering turned on. Note that including or not including the "hole" has 
only a small effect on the HBT results presented later. 

Figure [2] shows sample two-pion correlation functions from the rescattering model for r = 
0.1, 0.5, and 1.0 fm/c with fits to the Gaussian function, Eq. (5), for = 1.8 TeV p — p 
collisions. A comparison is also made for the r = 0.1 fm/c case between two px cuts on the 
pions, i.e. 0.2 < p T < 0.5 GeV/c andp^ > 1 GeV/c. As seen, the Gaussian fits qualitatively 
reproduce the trends of the model correlation functions, but do not well represent all of the 
details of the shapes, which include an exponential-like shape for r = 0.1 fm/c and some 
oscillatory behavior for r = 0.5 and 1.0 fm/c. The oscillatory behavior is a feature of the 
delta-function assumption of r, which becomes more prominent for larger values of r 

Figure [3] shows sample correlation functions where the model is run with a uniform 
distribution of r as a test, for the two cases r < 0.2 fm/c and r < 1.0 fm/c. As seen by 
comparing Figures [2] and El these cases closely resemble the correlation functions for r = 0.1 



8 



1(T 



i 

o 

c 

i 

f 



1000 



100 



PYTHIA p+p, E =1.8 TeV 



all particles, all p t 




hole — exclude 5% of n in 

-1 <y<1 and p >0.5 GeV/c 



a PYTHIA with hole 

o PYTHIA with hole+rescatt 

• PYTHIA 



§ 



-4 




*1 



1000 



c 

Is 10 



T3 



0.1 



PYTHIA p+p, E^- 1 .8 TeV 
all particles, -1<y<1 

hole — exclude 5% of n in 

-l<y<1 and p >0.5 GeV/c 



e 

A g 

A 



A PYTHIA with hole 

O PYTHIA with hole+rescatt 

• PYTHIA 



1 



P T (GeV/c) 



FIG. 1: Pseudorapidity and pt distributions from PYTHIA for p — p collisions at y/s = 1.8 TeV 
for three cases: 1) directly from PYTHIA, 2) PYTHIA with f] — pT "hole", and 3) same as 2) but 
with rescattering turned on and r = 0.1 fm/c. 



and r = 0.5 fm/c, respectively. A more complete comparison is shown later. 

Figure HH6] show comparisons of Gaussian fit parameters for pions between Tevatron data 
(Experiment E735) and model predictions with and without rescattering at r = 0.1,0.5 
and 1.0 fm/c, respectively, versus px and Nq. As anticipated earlier, rescattering is seen 
to have the greatest effect on the fit parameters for the smallest value of r, r = 0.1 fm/c, 
becoming less important as r increases until it is seen to have an almost negligible effect at 
t = 1.0 fm/c. All three r cases (with rescattering) do an adequate job of describing the flat 
dependence of A on px and Nc seen in E735. It is also seen that the overall trends of the 
E735 R dependencies, i.e. decreasing with increasing pt and increasing with increasing Nc, 
are best reproduced by the r = 0.1 fm/c case with rescattering turned on, the larger r model 
predictions becoming progressively flatter with increasing r. Another feature found in the 
correlation functions of the higher r values not found in the E735 correlation functions is 
the oscillation in the baseline seen in Figure [2J No such oscillation appears for the r = 0.1 
fm/c case, in agreement with E735. 

Figures [7] and [8] show results for running the model with flat distributions of r, i.e. r < 0.2 
fm/c and r < 1.0 fm/c, and comparing with the delta-function cases r = 0.1 and 0.5 fm/c, 
which are the average values of the two flat ranges, respectively, with rescattering turned 
on, and compared with E735. As seen, the fit parameters for the flat r distributions give 
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FIG. 2: Sample two-pion correlation functions from the rescattering model for r = 0.1, 0.5, and 
1.0 fm/c with fits to the Gaussian function for y/s = 1.8 TeV p — p collisions. 



virtually the same results as the delta-function r distributions, demonstrating that either 
method of running the model gives almost identical results. 

Summarizing this section, the main result of the comparison of the HBT fit parameters 
from the present model with those from Tevatron experiment E735 shown above is that the 
r = 0.1 fm/c case best describes the trends of the fit parameters on p T and N c . This would 
seem to imply that the hadronization time in these collisions is short, i.e. r << 1 fm/c, 
and, as a consequence of this, large hadron densities exist at the early stage of the collision 
resulting in significatnt hadronic rescattering effects. 
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FIG. 4: Comparisons of Gaussian fit parameters for pions between Tevatron data (Experiment 
E735) and model predictions with and without rescattering at r = 0.1 fm/c versus pr and Nq- 



B. Predictions for p-p collisions at = 14 TeV 

Figures M - EH present results from the present model calculations for p — p collisions at 
y/s = 14 TeV to simulate full-energy LHC collisions. All calculations shown are for pions at 
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FIG. 5: Comparisons of Gaussian fit parameters for pions between Tevatron data (Experiment 
E735) and model predictions with and without rescattering at r = 0.5 fm/c versus pt and Nc- 



mid-rapidity, i.e. — 1 < y < 1. Several values of r and cuts on particle multiplicity and pion 
Pt are studied. In the results shown below, multiplicity is defined as the total multiplicity 
of pions, kaons, and nucleons of all charge states for all rapidity and pt in a collision. In 
order to compare with future experiments, the approximate correspondence between the 
total multiplicity bins used and the more experiment-friendly average detectable (p? > 100 
MeV/c) mid-rapidity (—1 < y < 1) charged particle multiplicity is shown in Table I. Also 
shown in Table I is the fraction of minimum bias events corresponding to each multiplicity 
bin. From this it is seen that all multiplicities used are predicted to be easily experimentally 
accessible. 

The following choices were made for the conditions used in the model calculations in 
generating the LHC predictions: 

• Based on the comparisons presented above between the model and E735 results, pre- 
dictions using the "delta-function" model for r for the cases r = 0.1 and 0.5 fm/c were 
made. As shown above for the Tevatron calculations, these cases give almost identical 
results as for the "flat-distribution" model for r < 0.2 and 1.0 fm/c. Although the 
closest agreement between the present rescattering calculations and E735 results was 
obtained for the r = 0.1 fm/c case, predictions are also included for r = 0.5 fm/c 
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TABLE I: Approximate correspondence between the total multiplicity bins used in the present 
calculations and the average detectable mid-rapidity charged multiplicity. The fraction of minimum 
bias events is also shown. 
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FIG. 9: Pseudorapidity and px distributions from PYTHIA for p — p collisions at y/s = 14 TeV for 
two cases: 1) directly from PYTHIA, and 2) PYTHIA with r/ — pr "hole" and with rescattering 
turned on for the case r = 0.1 fm/c. 



since the hadronization time for y/s = 14 TeV collisions may be larger than that for 
y/s = 1.8 TeV. 

• Employ the same u 7]— hole" method as for the Tevatron calculations, i.e. the "hole" 
is defined by randomly throwing away 5% of the particles in the input PYTHIA events 
in the y — px region — 1 < y < 1 and pr > 0.5 GeV/c for the r = 0.1 fm/c case, and 
scaled down accordingly for the r = 0.5 fm/c case. Comparisons of pseudorapidity and 
Pt distributions between PYTHIA run for the maximum LHC energy and PYTHIA 
with the "hole" and rescattering turned on for r = 0.1 are shown in Figure [9j 

• To extract the pion HBT fit parameters from the invariant correlation functions gen- 
erated by the model, use Eq.(6). As described above, this is to better characterize the 
finer features of the correlation functions and thus get better fits to the calculations. 
Figures [10] and [11] show fits to sample model-generated correlation functions for the 
cases r = 0.1 and 0.5 fm/c, respectively. As seen, the fits are in general quite good. 

Figures [12] and Q~3] show the predicted dependences of the fit parameters on px and total 
multiplicity, m, for the cases r = 0.1 and 0.5 fm/c, respectively for y/s = 14 TeV p — p 
collisions. Plots are made with low and high multiplicity cuts, i.e. m < 100 and m > 300, 
and low and high p T cuts, i.e. 0.1 < px < 0.3 GeV/c and 0.9 < pt < 11 GeV/c. The 
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FIG. 10: Sample two-pion correlation functions from the rescattering model for r = 0.1 fm/c with 
fits to the general function (i.e. Eq.(6)) for ^/s = 14 TeV p — p collisions. 

behaviors of the fit parameters seen in Figures [12] and [13] are discussed separately below. 
1. R-parameter 



The .R-parameter, which is related to the "size" of the pion-emitting source (see Eq.(7)), 
is seen to have the largest variation for different kinematical cuts, i.e. the strongest depen- 
dencies on pt and m. These are greatest for r = 0.1 fm/c since they are almost completely 
due to rescattering effects. In Figure d2] it is seen that for the proper kinematical cuts R 
can be made to increase by a factor of three for increasing m or decrease by a factor of 
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FIG. 11: Sample two-pion correlation functions from the rescattering model for r = 0.5 fm/c with 
fits to the general function (i.e. Eq. (6)) for ^fs = 14 TeV p — p collisions. 

three for increasing p?. Experimental observation of such strong variations in R would be a 
convincing signature for the presence of significant rescattering effects and therefore a short 
hadronization time. 

2. X-parameter 

The A-parameter, which is related to the "strength" of the HBT effect, is seen to have 
weak dependences on the kinematical cuts and to have a similar magnitude for both the 
r = 0.1 and 0.5 fm/c cases. It tends to have a magnitude such that A < 0.5, which is 
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mostly due to the presence of long-lived resonances in the model calculations which tend 
to suppress it. Though weak, it has a slightly increasing tendency for increasing p? and 
decreasing tendency for increasing m, which is opposite the directions of the dependencies 
for R. 

3. a-parameter 

The a-parameter, which is seen in Eq.(6) to be related to oscillatory behavior of the 
correlation function, is seen to mostly depend on the value of r. As also seen in Figures [TU1 
and [TT| it is small, i.e. a < 0.2, for r = 0.1 fm/c and large, i.e. a ~ 0.5 — 0.6, for r = 0.5 
fm/c. The connection between a and r for the simple case of a delta- function hadronization 
time can indeed be shown to be r ~ a [9]. Experimental observation of oscillations in 
the correlation function, i.e. large a values, would be evidence for a larger value of the 
hadronization time, i.e. r > 0.5 fm/c. 

4- 5 -parameter 

As seen in Eq.(6), the ^-parameter is related to how "exponential-like" (8 ~ 1) or 
"Gaussian-like" (5 ~ 2) the correlation function appears. As with A, it is seen to have 
weak dependences on the kinematical cuts, to have a similar magnitude for both the r = 0.1 
and 0.5 fm/c cases, and to have dependencies which are opposite to the directions of the 
dependencies for R. It tends to have values in the range 0.7 < 5 < 1.5. 

IV. CONCLUSIONS 

A simple model assuming a uniform hadronization proper time and including final-state 
hadronic rescattering has been used to predict two-pion HBT fit parameters for p—p collisions 
at ^fs = 1.8 TeV and p — p collisions at 14 TeV. For small values of r, i.e. r < 0.5 fm/c, it 
is found that rescattering has a significant influence on the fit parameters. Comparing the 
model predictions with p — p experimental results at = 1.8 TeV, the closest agreement is 
found for small r, i.e. r ~ 0.1 fm/c. This suggests that 1) final-state hadronic rescattering 
is already important at -Js = 1.8 TeV and 2) hadronization times are short. 
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FIG. 12: General function fit parameters versus pt and particle multiplicity from the rescattering 
model with r = 0.1 fm/c for several multiplicity and px cuts for p — p collisions at y/s = 14 TeV. 
The dashed lines are drawn to guide the eye. 



As is seen in the above figures, there are significant differences in the magnitudes and 
dependences on kinematical variables for the general fit parameters evaluated at different r- 
values in the model calculations carried out for yfs = 14 TeV p — p collisions. Comparisons 
of these results with actual future data from the LHC will be able to establish a) if this 
simple model describes the data in even a qualitative way and b) if so, the scale of the 
hadronization time in these collisions. 
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FIG. 13: General function fit parameters versus px and particle multiplicity from the rescattering 
model with r = 0.5 fm/c for several multiplicity and pr cuts for p — p collisions at y/s = 14 TeV. 
The dashed lines are drawn to guide the eye. 
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